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<n: We consider cosmological scenarios based on f(R, T) theories of gravity (R is the 

Ricci scalar and T is the trace of the energy-momentum tensor) and numerically 
reconstruct the function f(R,T) which is able to reproduce the same expansion 
history generated, in the standard General Relativity theory, by dark matter 
and holographic dark energy. We consider two special f(R,T) models: in the 
first instance, we investigate the modification R + 2f(T), i.e. the usual Einstein- 

U 



Hilbert term plus a f(T) correction. In the second instance, we consider a 
f(R) + XT theory, i.e. a T correction to the renown f(R) theory of gravity. 

PACS numbers: 98.80.-k, 04.50.Kd, 95.36.+X 
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I. INTRODUCTION 

The nature of dark matter and dark energy is one of the most important issues today in physics. 
^ ■ There are strong observational evidences in astrophysics and cosmology for the existence of these 

two components of the cosmic energy budget, indicating that about 95% of the Universe is composed 
by dark matter (about 25%) and by dark energy (about 70%), but no direct detection has been 
reported until now. The usual candidates for dark matter (neutralinos and axions, for example) 
and dark energy (cosmological constant, quintessence, etc.) lead to very robust scenarios, but at 
same time they must face theoretical and observational issues. For recent reviews on the subject, 



see for example 



id- 



The strongest issue is perhaps the one regarding dark energy as the vacuum expectation value 
of some quantum field, which would be a natural candidate, but whose correct theoretical value 
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could be predicted only in the framework of a complete theory of quantum gravity, which still we 
do not possess. Nevertheless, it is possible, at least, to guess some of the features of this theory. 



In particular, the holographic principle 
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1 3f | may shed some light on the dark energy problem. 
According to this principle, in presence of gravity the number of the degrees of freedom of a local 
quantum system would be related to the area of its boundary, rather than to the volume of the 
system (as expected when gravity is absent). Following this idea, in [ul ] the authors suggested 
an entanglement relation between the infrared and ultraviolet cutoffs due to the limitation set by 
the formation of a black hole, which sets an upper bound for the vacuum energy. We can then 
interpret the ultraviolet cutoff as the vacuum density value, but still we need an ansatz for the 
infrared cutoff. As a candidate for such distance, in [l^i 1(J the authors propose and investigate 
the future event horizon, tested against type la supernovae data and cosmic microwave background 
anisotropies in 



17 



18( | . We shall present more detail on this in Sec. [UJ 



Adding new components of dark energy to the whole energy budget in order to explain the 
current observation is a way, but not the only one. Since General Relativity has been thoroughly 
tested up to solar system scales, it may be possible that the Einstein-Hilbert action contain correc- 
tions on larger, cosmological, scales thereby candidating as possible explanation of the evolution 
of the universe. Such modifications should be, in principle, composed by higher order curvature 
invariant terms (such as R 2 , R IL yR^ v ^ etc) but also by non-trivial coupling between matter or fields 



and geometry. See for example 



191423] for some reviews on the subject (especially on f(R) theory). 



It is also worth pointing out that these terms should naturally emerge as quantum corrections in 



26]. 



the low energy effective action of quantum gravity or string theory 

In this paper we connect these two approaches, considering a f(R,T) theory of gravity, where 
R is the Ricci scalar, whereas T is the trace of the stress-energy momentum. This modified 
gravity theory has been recently introduced in 27], where the authors derived the field equations 



and considered several cases, relevant in cosmology and astrophysics. As for the former, f(R,T) 
models have been constructed describing the transition from the matter dominated phase to the 
late times accelerated one [28]. 



Our task here, is to find out which form the function f(R,T) has to have in order to reproduce 



the same properties of the holographic dark energy proposed in 
the same reconstruction scheme proposed and employed in 



151 ] . To this purpose, we employ 



29j-[33|. For reference, in order to track 



the contribution of the T part of the action in the reconstruction, we consider two special f(R,T) 
models: in the first instance, we investigate the modification R + 2/(T), i.e. the usual Einstein- 
Hilbert term plus a f(T) correction. In the second instance we consider a f(R) + AT theory, i.e. a 
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T correction to the renown f(R) gravity. In both cases, we consider dark energy accompanied by 
a pressureless matter component (which would determine T). 

The paper is organised as follows. In Sec. HU the equations of motion are established and the 
holographic dark energy introduced. In Sec. IIIII and IIVI the above mentioned cases are analysed. 
Finally, Sec. [V] is devoted to discussion and conclusions. 

We use 8ttG = c = 1 units and adopt the metric formalism, i.e. the variation of the action 
is considered with respect to the metric quantities. 

II. f(R,T) GRAVITY AND HOLOGRAPHIC DARK ENERGY 

In [2?!, the following modification of Einstein's theory is proposed: 

S=\j f(R, T)^~g d A x + J L m ^g- d 4 x , (1) 

where f(R,T) is an arbitrary function of the Ricci scalar R and of the trace T of the energy- 
momentum tensor, defined as 

rp 2 S (^gL m ) 

^ sir ' 

where L m is the matter Lagrangian density. We assume the matter lagrangian to depend on the 
metric, so that 

Tftv = 9nvL m - 2 ^— . (3) 
Varying action ([1]) with respect to the metric g^ v , one obtains [3] 

f R (R,T)R^ - ±f(R,T) 9tu , + (g^B - V M V,) f R (R,T) = - f T (R,T)T^ - f T {R,T)Q^ , (4) 
where the subscripts RoxT imply derivation with respect that quantity and we have also defined 

Planning a cosmological application, we assume matter to be described by a perfect fluid energy- 
momentum tensor 

= (p + p) u^u u - pg^ v , (6) 

and that L m = —p, so that we have 

©aik = -2T^ -pg^ , (7) 
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and Eq. @ simplifies as 

/H(fl,T)12^-i/(12,T)^ + (^ . (8) 

In order to compare it with Einstein's, we cast the above equation as follows: 

f(R,T) 9flv -lg^R, (9) 



2f R (R,T) Jy ' 2 
where G^ u = R^ u — Rg^/2 is the Einstein tensor. Now we can identify 

as the effective matter energy-momentum tensor and 

T% eom) = - fR ( R T) " V M V„) T) + 2f J R T) f(R, T)g^ - \g, u R , (11) 

as the energy-momentum tensor of a "geometric" matter component. 

We now assume a background described by the Friedmann-Lemaitre- Robertson- Walker metric 

ds 2 = dt 2 - a(t) 2 5 tj dx i dx j , (12) 

with spatially flat hyper surf aces, and find a form for the function f(R,T) which is able to recon- 
struct holographic dark energy. 



A. Holographic Dark Energy 



According to the holographic principle 



11 



131 ] an entanglement relation between the infrared 



(IR) and ultraviolet (UV) cut-offs of a quantum theory, due to the limitation set by the formation 
of a black hole, sets an upper bound for the vacuum energy 



Pv 



36^ 
L 2 ' 



(13) 



where b is a free parameter and the IR (large scales) cutoff L needs to be specified by an ansatz 
We are interested in the one proposed in [151 ] : 

dt' r°° da 



L = R h 



(14) 



H a(t') J a H{d)d 2 ' 
i.e. the future event horizon, that is the distance covered by a photon from now until the remote 
future. Note that the very presence of a vacuum energy component makes the above integration 
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finite. We consider a model composed by holographic dark energy plus ordinary pressureless matter, 
i.e. 



3H 2 = p v + p m = p v + p m0 (l + zf , 



(15) 



where H = a/a is the Hubble parameter and the dot denotes derivation with respect to the cosmic 
time. Introduce the critical energy density p CT := 3H 2 , we define 

b 2 



Pv 



Per RlH 2 ' 

Using Eqs. (fT3|) and (fl4"|) . it is easy to show that 

b 



(16) 



R\< 



1 . 



(17) 



The holographic dark energy density p v evolves according to the conservation law 



p v + 3Hp v (1 + w v ) = 



(18) 



because in Eq. (|15p we have implicitly assumed the matter component to conserve separately. Now, 
using Eqs. (fT3|) . p^]) and (fT7|) . one can find 



b 



Rh Vv^v 



(19) 



Comparing (fT9]l with (|T8|) one can read off 



i 270; 

3 + 36 



(20) 



Moreover, combining Eq. (|M[) with Eqs. (|T3|) and (|T5j) . the evolution for J7 V is determined by the 
following equation: 



i 



i 



1 + 2 



ftv (1 - Hv) 



(21) 



where the prime denotes derivation with respect to the redshift z. Testing this model against 
type la supernovae and cosmic microwave background anisotropies, b turns out to be constrained 
around unity, with the case b < 1 favoured [l7l Jl^|. Note that, from Eq. (|20p . b < 1 means that 
the universe will end up in a phantom phase. For more comprehensive analysis of holographic dark 
energy models, we refer the reader to 



34, 



35 . 



In the next section we investigate a simple case of reconstruction of f(R,T). 
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III. A SIMPLE CASE 

We now consider a single perfect fluid model with density p and pressure p, together with the 
following ansatz (one of the first considered in [^tJ ) : 

f(R,T) = R + 2f(T) , (22) 

i.e. the action is given by the same Einstein-Hilbert one plus a function of T. This is a particularly 
interesting choice since, from Eqs. (jlOH and (llip . we get 

Tfj, u = (1 + 2/ T ) V + 2p/ T ^ + /(T)^ . (23) 

For p = one has T = p and, choosing /(T) = AT one can construct a model with an effective 
cosmological constant [3f3]. From Eq. ()23[) one can read off the effective energy density and pressure 
of the universe content: 

ZH 2 = p cS = (1 + 2/ r ) p + 2p/ r + f(T) , (24) 
-2tf - 3H 2 = PeS =p- f(T) , (25) 

and therefore a dark energy component may appear, even if we are considering a single perfect 
fluid model. From Eqs. (|24|) and ()25[) . it is clear that we can pick out a "fictitious" component, 
due to f(T), described by 

Pf = 2f TP + 2pf T + f(T) , (26) 
P/ = ~f(T) , (27) 

and, provided / positive, it may well describe a dark energy component, since its pressure is 
negative. In order to reconstruct the function / starting from the holographic principle, we note 
that the equation of state parameter of the dark component induced by / is 

"'--s&rniW (28) 

and we identify it with w v , the one provided by the holographic dark energy in Eq. (I20p : 

2(p + p)f T + f(T) 3 + 36 • 1 > 

For the standard model given in Eq. (I15p . consider the fluid component to be pressureless matter, 
i.e. p = 0. We are left to solve the following system of equations: 

H dp JyH, b + 2^/tl; v ' 

f = -h{ l+2J f) aA1 - a ^ (31) 
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where we have used T = p, because we are considering pressureless matter. In order to solve the 
above system, we have to fix some initial conditions. Clearly, Q v (p = po) = 1 — £l m o, and we choose 
^mo = 0.3, accordingly with current cosmological observation. For the initial condition on /, from 
Eq. d2U) (with p = 0) we write 

[l + 2/ T (p )]O m o + ^ = l. (32) 

Evaluating Eq. (f3"0"j) today and combining it with Eq. ([3"2"j) we find the following algebraic equation 
determining the initial condition on /: 

M j £$jH =4flJn -- (33) 

As we expected, when f2 v o = 0, then f(po) = and Eq. ([30]) implies that / is identically vanishing. 




0.02 -, 
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FIG. 1: Left panel: evolution of f(p). Right panel: evolution of f2 v . The cases considered are b = 
0.6,0.8,1.0,1.2 (solid black, dashed red, dot-dashed blue and dotted green, respectively). We have chosen 
as initial conditions in f2 m o = 0.3 the values J7 v o = 1 — fi m ii = 0.7 and fo given by Eq. (|33t . Note that / 
and p are normalised to 3ffg. The vertical lines in the plots represent p = po, i.e. the present instant. 

In FIG. [U we plot the solution of the system of differential equations (I30p and (]3ip . Note that 
we normalise / and p to 3Hq. As expected from Eq. (f30j) . for large values of p, i.e. far in the 
past, / oc p because the dark energy component is subdominant. The actual difference among 
the various choices of b takes place at late times, for small values of p. Again from inspection of 
Eq. (I30p . we can see that for large values of b the linear evolution / oc p is again solution. That is 
why in the left panel of FIG. Q]the curve seems to "straighten up" for increasing b. 
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FIG. 2: Same as FIG. [TJ but with the redshift z as independent variable. 



A final remark about the future evolution. It is clear from Eq. (|30p . that when J7 V — > 1, in the 
remote future, the solution for / gets the asymptotic form 

Hp) « p^ 2 • (34) 

Therefore, we would have a future singularity for —2 < b < 1, as it appears for the relevant cases 
of FIG. [T]and FIG. [5J The special case b = 1 implies an asymptotically constant /. 

We now turn our discussion on a more general case, where the curvature R comes into the 
action as a function to be determined. 



IV. A MORE COMPLICATED CASE 



Now we turn our attention to the special case 

f(R,T) = f(R) + \T , (35) 

i.e. a T-linear correction to the class of f(R) theories. With the ansatz (|35|) the matter content 
(fTU|) is "corrected" as follows: 

f -=^ r - + 7^) A! *- (36) 
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whereas the geometry induced stress-energy tensor is 

f% eom) = -j^r (g^D - V,*,,) f R (R) + ^—.\f{R) + AT]^ " \^R ■ (37) 

Our aim is now to reconstruct the form of the f(R) which is able to reproduce the holographic 
dark energy paradigm. We again consider a pressureless perfect fluid with density p and again 
assume metric (|12p . From Eqs. (|36p and (|37p the effective density and pressure are the following: 



r2 P , 3A R , / nrrfR 



m = p * = f R + y- R p —2 + W R -* H f R > (38) 

-2H - 3H 2 = p cS = J- (f R + 2Hf R ) ^ ■ (39) 



f R \ ' 2/r 2/r 2 

From Eq. ([38j) it appears that the energy density of the perfect fluid is rescaled by a factor 1/ f R . 
Looking at Eq. (|15j) . we can extract a form for p v in the following way: 

p « = Tr- , ' + w/-2 + w b -™t r - m 



whereas the form of p v is already given in Eq. ()39p . since our fluid is pressureless. From Eqs. ([39 
and ([40 p we can write the following differential equation for f R : 



h - Hf R -\p + Pv (i + w Y )} f R = -p(l + A) . (41) 

Note, as a cross-check, that for p v = A = the above equation simplifies to 

fR ~ Hf R - pf R = - P , (42) 

which possesses the particular solution f R = 1, i.e. / = R + A, the original Einstein-Hilbert action 
plus an integration "cosmological" constant. We expect this solution to be the only one, otherwise 
there would exist an alternative f(R) theory which would behave exactly as general relativity. Let 
us speculate a bit more on this point. If p v = A = 0, i.e. for a pure Einstein-de Sitter universe, we 
have from Eq. (j!5p 

ff=|, P = Po|, (43) 

where to is the present cosmic time (i.e. the age of the universe). Considering the homogeneous 
part of Eq. (|42|) and looking for a solution of the form f R oc t n , we find: 

n(n - 1) - -n - p t 2 = , (44) 

which gives: 



5 , 25 9 

"1,2 = - ± y — + p t 2 , (45) 
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and the general solution of Eq. (|42p can be written as: 



f R = i + Ci t ni + c 2 1 



(46) 



Now, the initial conditions we adopt here are /i?(to) = 1 an d /fl(io) = 0- The reason is essentially 



29H33J]. However, 



not spoiling the agreement between general relativity and solar system tests, see 
we stress here that these two conditions also imply C\ = C2 = and therefore restore the general 
relativity limit /r = 1 when p v = A = 0. 

Changing the variable to the redshift and employing Eqs. (|15ft and (I2ip one can recast Eq. (|41j) 
in the following compact form: 



(i + *)7S + 



1 + 2 



7 - o v 1 + 



i-a 



1 2VIX 
3 + 36 



-3(1 + A)(l-Jlv) 



(47) 



where again the prime denotes derivation with respect to the redshift z. The curvature R can be 
easily found as 



R = -6(i? + 2H 2 



3# 2 O m0 (l + zf 



1-fi, 

The initial conditions /ij(to) = 1 and //?(to) = 0, translated to the redshift variable, are 

d 2 R 

1 

z=0 



dz 2 



2=0 



dz 2 



dz 



2=0 



1 + J2v 1 



dR 

dz 



(48) 



2=0 



(49) 
(50) 



Finally, the initial condition on / can be extracted by Eq. ([38]) . being that p v o = 3H$ — p m o- Thus, 
we have 



f( z = 0) = r( z = 0) + ml (1 - n m0 - hn m0 ^j 



(51) 



Note the correction due to the AT term. In FIG. [3] we plot the solution for /. Note that we 
normalise / and R to 3Hq. In the left panel we fix A = 0, i.e. we are actually considering a pure 
f(R) theory, and vary b. In the right panel, on the other hand, we consider positive and negative 
values of A. As one may note, A has a poor influence on the evolution of /. We could expect 
this from inspection of Eq. (|47p . Indeed, A only enters the source term on the right hand side and 
therefore, when S7 V grows to unity, its impact on the evolution of / is weak. On the other hand, 
larger values of A may have a relevant effect at early times, determining the slope of /. 
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FIG. 3: Evolution of / as a function of the curvature R. Left panel: A = 0, i.e. we are considering 
a pure f(R) theory, and b — 0.6,0.8,1.0,1.2 (solid black, dashed red, dot-dashed blue and dotted green, 
respectively). Right panel: b = 1.0 and A = —0.2,0,0.2,0.4 (solid black, dashed red, dot-dashed blue 
and dotted green, respectively. The curves appear superposed.). We have chosen as initial conditions in 
fl m o = 0.3 the values f2 v o = 1 — ^mo = 0.7. Note that / and R are normalised to 3Hq and the redshift 
interval chosen is < z < 10. 

In FIG. [5] and FIG. [5] we display the future evolution of /, both as function of R or of z. For 

the same reason stated above, the effect of A is not relevant. 
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FIG. 5: Same as FIG. IU but with / as function of z. 

Finally, in FIG. [6] we plot the solution for Q v . Again, its evolution appears to be independent 
of A. 
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FIG. 6: Evolution of Q v as a function of the redshift. Left panel: A = 0, i.e. we are considering a pure 
f(R) theory, and b = 0.6,0.8, 1.0, 1.2 (black, dashed red, dash-dotted blue and dotted green, respectively). 
Right panel: b — 1.0 and A = —0.2,0,0.2,0.4 (black, dashed red, dash-dotted blue and dotted green, 
respectively). Note that the curves are superposed. We have chosen as initial conditions in il m o = 0.3 the 
values ri v o = 1 — fl m o = 0.7. 

V. DISCUSSION AND CONCLUSIONS 

In this work we have investigated a description of holographic dark energy in terms of suitably 
reconstructed f(R,T) gravity theories. The latter have been recently introduced as modifications 
of Einstein's theory possessing some interesting solutions which are relevant in cosmology and 

n 

astrophysics [271 ] . 

We have considered two special types of models: f(R,T) = R + 2/(T), i.e. a correction 
to the Einstein-Hilbert action depending on the matter content, and f(R,T) = f(R) + AT, i.e. a 
simple T- linear correction to the class of f(R) theories. 

Since we have assumed the matter content to be a pressureless perfect fluid, then T = p, i.e. 
the corrections assumed are directly dependent on the energy density of the universe content. We 
have constructed differential equations for the function / under investigation and numerically 
solved them, physically specifying the required initial conditions. Our simple analysis shows that 
holographic dark energy models are contained in the larger class of f(R,T) theories, at least 
considering a given background evolution of the universe. 

It would be interesting to investigate how the evolution of matter perturbations would change, 
depending on the description of dark energy. We expect, in principle, different results when using 
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holographic dark energy or its f(R, T) reconstruction and therefore there is possibility for discrim- 
inating between the two descriptions. For example, it would be interesting to adapt the recently 
proposed scheme for perturbations in f(R) theories |37l] to the broader f(R,T) class. We leave 
this as a future work. 
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